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Abstract 

Being a determination at low energies, the analysis of hadronic t decay data provides a rather precise determination 
of the strong coupling a, after evolving the result to Mz- At such a level of precision, even small non-perturbative ef- 
fects become relevant for the central value and error. While those effects had been taken into account in the framework 
of the operator product expansion, contributions going beyond it, so-called duality violations, have previously been 
neglected. The following investigation fills this gap through a finite-energy sum rule analysis of t decay spectra from 
the OPAL experiment, including duality violations and performing a consistent fit of all appearing QCD parameters. 
The resulting values for q's(Mt) are 0.307(19) in fixed-order perturbation theory and 0.322(26) in contour-improved 
perturbation theory, which translates to the nj - 5 values 0.1169(25) and 0.1187(32) at Mz, respectively. 
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1. Introduction 

The presented contribution summarises the recent ar- 
ticle 1 1 1, which laid out a new framework to determine 
the strong coupling a, from hadronic r decays. After 
the determination of the order coefficient to the QCD 
Adler function |2|, interest in the extraction of a, from 
T decays was revived ||2l [3] S |5] . However, the results 
obtained in two approaches to resumming the perturba- 
tive series, fixed-order perturbation theory (FOPT) and 
contour-improved perturbation theory (CIPT), are only 
barely compatible. 

While this issue will be left untouched here, the 
determination from r decays also requires numerical in- 
puts for the non-perturbative contributions that arise in 
the framework of the operator product expansion (OPE), 
the QCD condensates. Even though the OPE terms are 



usually included in the analysis, contributions which go 
beyond it, so-called duality violations (DVs), and which 
arise due to the presence of resonances on the physical, 
Minkowskian axis, have so far been neglected. How- 
ever, for a consistent estimate of the OPE contributions, 
effects of DVs should be included since typically several 
spectral integrals of the experimental data with differing 
weight functions, and possibly at varying energies, are 
employed, to which the DVs contribute differently[^ An 
effort to fill this gap was initiated in the work of ref. HI. 

2. Theoretical framework 

The central observables for the analysis of hadronic t 
decay are the total decay rate into hadrons, 

_ r[T" ^ Vt hadrons] 
Rt = 1 ' U) 

T[t -^VrC Vj 
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and corresponding differential decay spectra. Experi- 
mentally, those spectra can be decomposed into the non- 
strange vector (V) and axial-vector (A), as well as the 
strange components. 

The theoretical description is based on the two-point 
correlation functions 

U'^'^^iq) = ij d'xe"'^0\T{jyx)r^JOm} (2) 

where the relevant V and A currents are given by 
Jy(x) - u(x)y^d(x) and J'^ix) - u(x)-Y''ysd(x), and the 
superscripts (1,0) label the spin. 

Employing the physical spectral functions py\{s) = 
ImTly^^is)/?:, the V and A contributions to eq. |lj can 
be expressed as the sq = version of the weighted 
integral Q 



Rv.a(so) = Utt^SewW, 



ud\ 



J So \ So' 



50 / So ' 



(3) 



The scalar/pseudoscalar (J - 0) terms are suppressed 
by factors of and hence negligible, apart from the 
pion-pole contribution. Thus, the experimental decay 
distributions essentially determine py^^\s). 

Making use of Cauchy's theorem, the analytic struc- 
ture of Y]^y^''\s) implies that for any analytic weight 
w{s) the following finite energy sum rule (FESR) is sat- 
isfied 18]: 



dswis)I\^^fHs).iA) 



l-s|=io 



The central idea then is to evaluate the LHS of Q from 
experimental data, and the RHS within the theoretical 
QCD description, thereby allowing to determine funda- 
mental QCD parameters like ff ,. 

Conventionally, the theoretical side is calculated in 
the framework of the operator product expansion. How- 
ever, the contour integral has to be performed down to 
the physical, Minkowskian axis, where the OPE breaks 
down due to the presence of resonances, and quark- 
hadron duality violations arise. Even though the kine- 
matic weight in (|3]l has a double zero at i = sq, which 
may suppress DVs, for a general weight wis) those con- 
tributions should be included. 

The theoretical 2-point correlators Xl^J'^Vi) will thus 
be represented as 



where the second term corresponds to the DV contri- 
bution. Further assuming that it vanishes sufficiently 
fast at infinity, the required contour integral containing 
Ay^(i) can be expressed as 



2n 



CO 

j)dsw{s)^vA{s) = - J dsw{s)p'^l(s). (6) 



|i|=io 



Guided by a model for the light-quark V/A correlators 
im, being based on Regge theory and large-A^c, we em- 
ploy forp^^(i) Ansatz IfTOllTTI 



Pvj^is) = KvA^ sin (ay,^ +I3v.as) ■ 



(7) 



U'^f\s) = nZ%) + Av,A(s), 



(5) 



Such a model for DVs is also supported by a study of 
the Coulomb system 1 12 1. Besides the QCD parameters, 
per channel this adds four additional parameters into the 
theoretical description, all of which should be extracted 
through fits to experiment. 

3. Numerical analysis 

The central strategy to extract all parameters, the 
strong coupling as, the OPE parameters, as well as DVs, 
is a simultaneous fit to weighted spectral integrals (mo- 
ments) of the form Q with in general different weights 
w(s) and/or different so, with several criteria restricting 
the set of sensible choices. 

First of all, for simplicity, only weights polynomial 
in s will be considered. Next, as we intend to ex- 
tract the parameters for DVs, a weight without pinch- 
suppression (a zero at 5 = so) should be included. The 
simplest choice here is w(s) = 1. The order of the se- 
lected polynomial also determines the dimension up to 
which OPE terms contribute at leading-order (neglect- 
ing logarithmically suppressed effects). For a poly- 
nomial of highest order s", operators up to dimension 
D - 2n + 2 contribute. As very little is known about 
OPE contributions beyond D - 8, w(i) will be re- 
stricted to at most 3rd order The kinematical weight 
of eq. ([3]) is of this type. Furthermore, below some min- 
imal So ~ 1.5 GeV^, both, the OPE as well as the per- 
turbative expansion become questionable and also the 
DV model is no longer adequate. Hence only moments 
down to such a value should be included. Finally, the 
moments for considered weights/so combinations are 
correlated, and for the fits to work these correlations 
should not become too strong. This is the most delicate 
point about our analysis. 

The required V and A spectral functions have been 
experimentally determined by the ALEPH [13, 14J and 
OPAL ifTSl collaborations at LEP. However, the most 
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recent publicly available ALEPH data lfT4l do not fully 
include correlations due to unfolding. Therefore, we 
chose to perform our analysis on the basis of the OPAL 
data IBj only. 
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Figure 1: Comparison of experimental and tlieoretical w(s) = I V 
moments. FOPT fits are sliown in blue (solid) and CIPT fits in red 
(dashed). The (much flatter) black curves display only the OPE parts 
of the FOPT (solid) and CIPT (dashed) fit results. 

Let us begin by discussing the most basic fit employ- 
ing only w(s) = 1 as the weight. In this case, sq's corre- 
sponding to the upper endpoints of all experimental bins 
down to imin ~ 1 .5 GeV^ can be included if the fit is per- 
formed only on the V spectrum, and the fit parameters 
then only consist of a, and the four V-channel DV pa- 
rameters. The comparison of the moments as computed 
from the OPAL data and theoretical prediction with fit- 
ted parameters is displayed in figure 1 . The blue (solid) 
and red (dashed) lines correspond to FO and CI pertur- 
bative results respectively. For comparison the flatter 
black solid and dashed lines correspond to the same fits 
omitting DVs, which clearly demonstrates that it is nec- 
essary to include them in order to obtain a reasonable 
fit. The a.5 values resulting from the fit read 

a,(MT) = 0.307 + 0.018 ±0.004 (FOPT), (8) 
as(Mr) = 0.322 + 0.025 ± 0.004 (CIPT) , (9) 

where the first error is the fit uncertainty and the second 
results from a variation of i^in. The corresponding DV 
parameters can be found in table 1 of [H^. Similarly to 
previous ofj determinations from r decays [2, 3^ 4] |5l, 
the FOPT result turns out to be lower than the CIPT 
value, though due to the additional DV parameters the 
errors turn out larger and the difference is less signifi- 
cant. 

The next level of sophistication is including the axial 
channel in the fit with w{s) = 1. This introduces the 
axial DV parameters as additional fit parameters. It is 
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Figure 2: Comparison of experimental and theoretical H'(i) = 1 A 
moments. FOPT fits are shown in blue (solid) and CIPT fits in red 
(dashed). The (much flatter) black curves display only the OPE parts 
of the FOPT (solid) and CIPT (dashed) fit results. 



found, however, that a fit employing the full set of sq, 
corresponding to all right bin endpoints above is 
no longer possible. A stable fit is still obtained, though, 
if only every third so is included. The comparison of 
the vector moments in this case looks very much like 
figure 1. In figure 2, the corresponding comparison is 
given for the A moments, with the notation being the 
same as for figure 1. Again, a fit without DVs would 
not provide an acceptable description of the data. The 
resulting values for a, are found to be 

as{M,) = 0.308 ± 0.016 + 0.009 (FOPT), (10) 
as{M,) = 0.325 + 0.022 ±0.011 (CIPT), (11) 

very similar to the vector only case, but somewhat less 
stable under variation of i^in. A full account of all fit 
parameters is found in table 3 of |[T|. 

Regarding fits including several weight functions, the 
additional weights w(x) - 1 - with x = s/sq and the 
kinematical weight w(x) = (1 - x)^(l + 2x) have been 
investigated by us in detail. The first one has a single 
pinch suppression and includes the D - 6 OPE term 
while the second is doubly pinched and involves con- 
densates up to D = 8. As soon as a second weight is 
added to the fit, even reducing the number of employed 
sq only leads to very unstable fits, due to the strong cor- 
relations. Thus, for these fits we have decided to fol- 
low a different route. Here the cross-correlations be- 
tween different weights are dropped in the fit, but later 
again included by a linear fluctuation analysis. This 
again yields results fully compatible with the ones given 
above, though without a reduction in the error of o-j. 

As one example figure 3 displays the moment com- 
parison for the kinematical weight. Due to the dou- 
ble pinch suppression, the pure OPE fares much bet- 



/ Nuclear Physics B Proceedings Supplement 00 (2011) 7-^ 



4 



0.019^ : 

0.018 ; 

0.016 -f ^ ^^^^'iSitoiiii. 
0.015- / 

0.014 i/ : 

o.ois j/ 

1.0 1.5 2.0 2.5 3.0 

so (GeV^) 

Figure 3: Comparison of experimental and theoretical w(x) = (1 - 
x)^(i + 2x) V moments. FOPT fits are shown in blue (solid) and CIPT 
fits in red (dashed). The black curves display only the OPE parts of 
the FOPT (sohd) and CIPT (dashed) fit results. 



ter even down to lower energies. Still, below roughly 
So ~ 2.5 GeV^ differences to the full fit including DVs 
are clearly visible. Furthermore, the inclusion of the 
DVs also influences the values of the fitted condensate 
parameters. This is exemplified by a more detailed dis- 
cussion of the D - 6 condensates and violations of the 
vacuum-saturation approximation which can be found 
in section VLB of ref. HI. 



Such an analysis is currently under way. Then, for com- 
parison, an analysis of the eventually updated ALEPH 
data would be most helpful, especially as the errors from 
this data set are expected to be smaller. Finally, in the 
long run, it is to be hoped that also the B-factories BaBar 
and Belle will at some point provide the complete V and 
A spectral functions, since with the available statistics, 
substantial improvements over the LEP experiments are 
to be envisaged. 
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4. Conclusions 

Based on ref. 1 1 i, a brief summary of the determi- 
nation of el's from hadronic t decays including effects 
due to duality violations in a consistent fashion has been 
presented above. The central result will be taken from 
eqs. ([8| and (9), since as far as OPE contribution is con- 
cerned, the weight w = 1 is most clean. Evolving this 
result to the Z boson mass, one obtains: 

a,(Mz) = 0.1169 ±0.0025 (FOPT), (12) 
a,(Mz) = 0.1187 ±0.0032 (CIPT). (13) 

Our results are somewhat lower than the most frequently 
quoted previous values from t decays, but on the 
other hand the uncertainty is increased in view of the 
additional degrees of freedom in the fit through the DVs. 

The analysis presented in fT] and this contribution 
was based on the original OPAL data [151, in order to 
clearly be able to compare to this analysis and single out 
the effect of including DVs. The corresponding results 
by OPAL for aAMz) read 0.1191(16) and 0.1219(20) 
for FOPT and CIPT respectively. 

As an obvious next step, the analysis should be re- 
peated with an updated data set which includes up-to- 
date T branching fractions and values of other inputs. 
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